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ABSTRACT 

We propose and analyze several block iteration preconditioners for the solution of elliptic prob- 
lems by spectral collocation methods in a region partitioned into several rectangles. It is shown that 
convergence is achieved with a rate which does not depend on the polynomial degree of the spectral 
solution. The iterative methods here presented can be effectively implemented on multiprocessor 
systems due to their high degree of parallelism. 
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INTRODUCTION 

We consider the numerical solution of elliptic problems by spectral collocation methods in regions 
partitioned into several rectangles (also called subdomains). The interaction between the 
subdomains is handled by proper iterative methods. 

Such iterative domain decomposition methods are particulary attractive in the context of spectral 
approximations. Actually, they allow the reduction of a problem set in a complicated geometry to a 
sequence of problems of similar type but with smaller size in every subdomain. The latter can 
therefore be faced by standard monodomain spectral solvers. Besides, due to their (hopefully high) 
degree of parallelism, these methods can be advantageously implemented on multiprocessor 
systems. 

To be effective, the rate of convergence of the iterative method should not deteriorate as the 
polynomial degree N of the numerical solution in each subdomain increases. This is what we prove 
for our method which alternates the solution of Dirichlet problems on the odd subdomains to that of 
mixed Neumann-Dirichlet problems on the even ones. Our proof applies to any domain Q 
decomposed into rectangles without internal crosspoints, i.e., whose vertices lie on the boundary 
of Q (see ftg.0.1). In the last section we consider a decomposition with internal cross points; in 
this case it is possible to prove the convergence of the iterative method for fixed N, although the 
rate of convergence decreases if N increases. This convergence behavior with respect to the 
decomposition is shared by other multidomain approaches based on finite elements or finite 
differences. We refer, e.g., to [4], [5]. 
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FIG.0.1 a) Decomposition without internal cross points b) Decomposition with an internal cross 
point 

The results we obtain generalize those of [15] which were established with a different proof 
technique and were limited to the case of a rectangle divided into two subdomains only. 

The outline of this paper is the following. In section 1 we write the differential problem in a 
variational form which is set up on the multidomain decomposition of the physical domain. 
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In section 2 we introduce the multidomain collocation problem. This amounts to collo- 
cating the differential problem at the Chebyshev (or Legendre) collocation points internal 
to each subdomain, and to enforcing, at the interface-boundary points, the continuity of 
the solution as well sis of a suitable combination of the residual of the equation and the nor- 
mal flux. This method has been proposed by Funaro [13] and retains the same convergence 
properties of the more classical patching method introduced by Orszag [23]. In particu- 
lar, as shown in [13] and [19], the numerical solution converges with spectral accuracy to 
the physical solution as the polynomial degree N grows on each subdomain. In section 2 
we also define an “iteration-by-subdomains” method for solving the discrete multidomain 
decomposition problem. Each step of this iterative method is based on the solution of 
Dirichlet problems in the odd subdomains, and mixed problems in the even ones. In sec- 
tion 3 we give a variational formulation of both the multidomain decomposition problem 
and the iteration-by-subdomain method. These variational formulations allow us to prove, 
in section 4, the convergence of the iterative method, the rate of convergence being inde- 
pendent of the degree of the polynomial solution in each subdomain. In section 5 we show 
that the iteration-by-subdomain method here analyzed amounts to solving the capacitance 
system governing the interface unknowns by a linear iteration procedure using a proper 
block diagonal preconditioner. Going further along this equivalence, we propose in section 
6 other preconditioners for the capacitance system, which in turn gives rise to some new 
iteration-by-subdomains algorithms. Finally in section 7 we consider a decomposition of 
the domain with internal cross points and we state how our method can be formulated 
and analyzed in this case. 

Numerical experiments on this iteration-by-subdomain method are presented in [13], [28]. 

Earlier works on algorithms of this nature for finite element and finite difference domain 
decomposition methods are reviewed and analyzed in [26]. 
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1. MULTIDOMAIN FORMULATION OF THE DIFFERENTIAL PROBLEM 


Let ft be an open two dimensional domain with boundary 3 £2. We consider the 
boundary-value-problem 

(1.1) Lu = f in £2 , u = 0 on 8£2 
where f is a given function of L 2 (£2) and 

(1.2) Lu:=-Au + a 0 u , a Q (x) £ 0 Vxe£2 , a Q e C°(£2) 


In variational form (1.1) reads as 


(1.3) 


u€Hj(Q): a(u,v): = 


(Vu Vv + a 0 uv)dx = 


fvdx V v e H{, (£2) 


n 




The Lax-Milgram theorem ensures that there exists a unique solution to (1.3); moreover, 
ue H 2 (£2), hence (1.1) holds almost everywhere in £2 (see, e.g., [18]). 

Let us consider a partition of £2 by non intersecting open subdomains £2 if i = 1,...,M, as in fig. 1.1, 

and denote by F the common boundary between £2 ; and £2 i+1 , i < M - 1. Then set T = UT. and, 

M-l 

for each i, we denote with <b := n (T.) the space of functions defined on T, which are the 

i=l 


traces on T of functions of H 1 0 (£2) (see [18]). Note that these functions vanish at the endpoints of 
Tj for i=l,...,M-l. 



Fig. 1.1: The partition of £2 
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For any <peO, let <j> e H ! o(0) denote any extension of 9 to Q, An example is furnished by the 
"harmonic" extension R<p, whose restriction Rj<p to Q, satisfies the Dirichlet problem: 


(1.4) 


aj(Ri<p,v) = 0 
^9 = 0 
R.<p = <p 


VveH^Qi) , 
on n 3Q , 
onTjUT. , 


where we have set 


ai (u,v) 


» 

(Vu Vv + ctjj uv) dx . 

J Qi 


The following equivalence statement is readily proven. 

Proposition 1.1 Problem (1.3) is equivalent to find ue Hyii) 59 , .that u ; := u^ , 1 S i < M, gr£ 
50lHti.QnS.iQ- 


(1.5) 


aj (Uj.v) = 


fvdx 




Vve Hj (Ity . 1 £i£M , 


(L6) Ui = u i+ i 

(1.7) Uj*0 


on T. , 1 £ i < M-l , 

on dSl^r>dQ , lSi^M , 


M 

. 


M 


/* 

(1.8) X • 

i*l 

a i( u i>9) * 1 
J 

ftpdx 

►*- X < 

j*2 

a/Uj.9)- 1 

fcpdx > 

i odd 



j even 

J 

n j 


v«p€<b . 


Qearly, (1.5) amounts to require that 

(1.9) LUj*f in 12; , l£i£M 


while (1.8) is equivalent to the transmission condition 


9uj 3u i+1 

(1.10) — = — onT , liliSJM-l 

do do 


(o is the outward unit vector to I". ). 
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2. SPECTRAL COLLOCATION APPROXIMATION TO THE MULTIDOMAIN PROBLEM 

In this section we consider a spectral collocation approximation of problem (1.5)-(1.8) 
which consists of collocating the differential equations (1.9) at Gaussian collocation points internal 

to Qj and to enforcing the interface conditions at some selected points of T. Either 
Gaussian-Legendre and Gaussian-Chebyshev points will be considered. We also give a variational 
formulation of the multidomain problem, then we define an "iteration-by-subdomains" procedure to 
solve effectively this problem. 

We confine ourselves to the case of a rectangular domain Q partitioned into M non intersecting 
rectangles without internal cross points (see fig.2.1), but what we are proving can be extended 
to more general situations as well (cf. remarks 3.1 and 4.2). 
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FIG.2.1 Decomposition by aligned subdomains 


Let D be the reference domain (-1,1) 2 and let N x , Ny be two given natural numbers. The Legendre 

Gauss Lobatto collocation points in D are the roots (£ km , 0$k< N x , 0 < m < Ny) of the 
polynomial 

a a 
gj • 

where L^t) is the Legendre polynomial of degree k in [-1,1]. In the Chebyshev case, the Gauss 
Lobatto collocation points in D are given by 

„ ( Jtk jtm'v 

= [* COS FT*. - cos n“ J ’ for0 ^ k ^ N x . 0£m£N y . 
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Wc remind that, in both Legendre and Chebyshev cases, 2^ lies on the boundary of D if either k 

equals 0 or N x , or m equals 0 or Ny. Moreover 

C i : = (04.yJ.OaSkSNi.OSmS N y } , for i=l,...,M , 

will denote the set of collocation points in the domain Qj. Given M natural numbers N' x , i = 
we take the points of C s as the images of the Legendre (or Chebyshev) points in D, with 
N x = N l x , through the affine transformation which maps D into For the convenience of 
exposition we set 

C““ :« q n Qj , C* := q n dCl , C ' :» qycf' u cf) for i=l M . 

See fig.2.2 for an example with M=2. 

0 points of C^i 
n points of C b 2 

♦ points of C^ nt 

• points of C2‘ nt 
© points of C^ 1 
0 points of C^j n C b 2 

Ftg.2.2 Interior and Boundary Collocation Points 

For i=l,...,M we denote by P N (Clj) the space of polynomials of degree N ! x with respect to the x 
variable and degree Ny with respect to the y one. Moreover we set 

P° (Qj):-(pe P N (£2 i ):p = 0onaa i ) . 



For the convenience of the reader we recall the Legendre and Chebyshev Gauss-Lobatto formulae 
that will be extensively used in the sequel (see also [10]). 


( 2 . 1 ) 


•7 - 

Lesendre Gauss-Lob atto formula 

*4 N y 

g(x,y)dxdy=£ L , y m ) co‘ k co m , 

k=0 m=0 

Oi 


where 


( 2 . 2 ) 


<v 


a i 


i m2 


Ni(N* + l) [M<)] 


and co = 


m N y (N y + 1) 




We have set aj * x { - x iA , where Xj denotes the abscissa of H. 

Chebvshev Gauss-Lobano formula 


(2.3) 


K Ny 

g(x,y)co (x,y)dxdy=S Z g(*k.y m ) co' k co m , 


k=*0 m*0 


where 


(2.4) 


0) (x,y) * [(x - x i . 1 )(x i - x)]’ 1/2 (1 - y 2 y m , 

i a*7C * i 

ca * — r for k=l,...,N - 1 , co. = — r for k=0 and k=N‘ , 

K 2JSJ 1 K vIVT^ * 


4N: 


Jt it 

C 0 m = rr- for m=l,...,N v - 1 , co = -stt- for m=0 and m=N„ 
m N y y m 2N y y 


We can now introduce the spectral collocation approximation to the problem (1.1), which 
fits well the multidomain formulation (1.5M1.8). The spectral solution u N verifies: 


(2.5) 

( 2 . 6 ) 

(2.7) 

(2.8) 


U IN u NI£lj 6 P N^i) » i = l.--,M , 


Lu^ = f 

inCf‘ , i=l,...,M , 

U IN ~ U i+1,N 

n 

inC , 

U i.N “0 

inC|* , i=l M , 

^ U i+1,N 

n — ^r = 

- (Lu^ - 0 0)j - (Lu i+ , 


n 


(2.9) 
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* i + i+1 

Here we have set to. = co . and to. = o> 0 for We are assuming f to be a continuous 

function in order that (2.6) and (2.9) make sense. Note that at the interface collocation points, the 
jump of the normal derivative is not asked to vanish (as in (1.10)), but rather to balance a suitable 
linear combination of the residual of the equation from both sides. We observe,however, that the 
coefficients of such combination tend to zero if N* x tends to infinity. This follows from the 

definitions of the weights co' k of the Gauss-Lobatto formula given in (2.2) and (2.4). 

From a computational point of view, dealing with condition (2.9) rather than with the pure flux 
d d 

condition u^ = g^ VijJ does not any appreciable extra work. On the other hand, a 

considerable advantage is achieved for the theoretical properties of the scheme, as is 
will be pointed out in the forthcotr 


Remark 2.1 The collocation scheme (2.5M2.9) has been first proposed by Funaro [13], who 
stressed its equivalence with a suitable variational approximation to (1.5)-(1.8) (see below).This 

equivalence allows one primarily to investigate the convergence properties, for N‘ x , N y -» <», of 
the method. 

The most relevant result (for the Legendre collocation) is the following (see [13] and [19]). 

0 . 

Setting Nj = min(N x , N y ) and assuming that Uj € H ‘(Qj) with c. £ 1, for i=l,...,M, the following 
error estimate holds 


( 2 . 10 ) 


M 

Uj * Ujjtf II , N:" 0 -* (II U: II + II f II , ) 

1 w* J Jo- o-l 

H (Qi) J“ l H J (tip H 1 (tlj) 


where C is a positive constant independent of Nj and uj, for j=l,...,M. In particular this estimate 

shows that if f = 0 then u i N = 0 for i = 1, .... M, hence the discrete problem has a unique 
solution. 

Another implication of (2.10) is that as N* x , N «o the jumps in the normal derivatives 

a a 

5x U »- N * 5x u i+i.N at t " e interface points vanish with spectral rate, hence so does the right hand side 

of (2.9). In [13] it is also shown by numerical experiments that using (2.9) rather than the pure 
flux condition does not reduce the accuracy of the results. 0 

We now introduce an iteration-by-subdomains procedure for the solution of the above 
collocation problem. Let us suppose, here and in the sequel, that M is an odd number, the case of 

M even can be studied analogously. We set T = u T , and define 
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P°N (T) := {p : r -> R : Pj := p| n is a polynomial of degree £ N y on I*., 
vanishing at the endpoints of F). 

For a given g°€ (T) we look for a sequence u^.n £ 1 , satisfying: u^e i= 1 . .,M,and: 

faduodd: 


(2.11) Lu"^ = f 


inCj®' , 


(2.12) uj'jj = g" 4 in C* 4 u C* for i * 1,M, u" 1N = g nl in C\ uJ, N = g nl in C M \ 


(2.13) u'Vn-O in C b j ; 

fori even: 

(2.14) Lu 11 ^ = f in 


(2.15) 

(2.16) 


*Em 


+ ^ U M.N 


Kn 




+ (L U ^-f)<O i »-g^-(L< 1 -f)<D i 



in C 

+ 1 

f) CO. 

in C 


i-1 


(2.17) u\n = 0 inC b i 




Then set: 




0 u" - N + (1 - 9 ) g" 4 

(2.18) g" = - " * 1,N 

on C ' A 

for i=3,...,M 

, i odd 


T - 

onC 1 

fori=l,...,M-2 

, iodd 


and restart from (2.1 1) with n+1 instead of n. In (2.18), 0 n is a positive relaxation parameter, how 
to choose it will be dealt with in the remark 5.4. 

Following this method, at each step one solves the M/2+1 independent Dirichlet problems 
(2.11)-(2.13); then, after computing the (independent) right-hand sides of (2.15)-(2.16), one has 
to solve the M/2 independent mixed problems (2.14)-(2.17). The degree of parallelism of this 
algorithm is therefore M/2. 

In the simple situation of a decomposition bv two subdomains , the above method simplifies 
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to: 


(2.19) 

Lu"i = t 

(2.20) 

“% “g"' 1 

(2.21) 

=° 

(2.22) 

(2.23) 

LuHjjj = f 
du" N 

-^--(Lu^-f)co 

(2.24) 

u\n = 0 


in Cj” 11 , 

inCj r 

inC 11 ! , 
inCj™ 

9uj N . _ 

- jj x + (Lujjj - f) cOj in C 

in C b 2 


and g" = 0 „u n 2 jM + (1 - 0 n ) g"' 1 on C r . Here T is the interface between Qj and 


Remark 2.2 (More general differential operators! The multidomain approach (2.5)-(2.9), and the 
corresponding iteration-by-subdomain procedure (2.11)-(2.18), apply also to the more general 
differential operator 


(2.25) 




+ a 0 u 


with a.j symmetric, uniformly negative definite, bounded, and smooth on £2. In such a case, the 
operator L in Qj should be replaced by its pseudo-spectral realisation 


(2.26) 


r T i f 3u 

3u Y 

1 d l\i 

i ( du 

3uV 


LH°n3£ 

+ a i2^J 

J + 5y L I] 

NV a 21^ + Ct 22 

Ty). 

+ a 0 u 


Here r N v denotes the polynomial of P N (£2j) which interpolates v at the collocation nodes of C ; (see 

[6], Ch. 2). Moreover, the conormal derivative a n 3/9x + a 12 9/3y should be used instead of the 
normal derivative 9/dx at the subdomain interfaces. 

In particular, we remind that this situation occurs whenever we consider the Laplace equation 
within a domain £2 subdivided into subdomains with curved boundaries. The mapping of the 

physical subdomain Qjinto the computational one D=(-l,l) 2 introduces a transformed operator 
which has the form (2.25).0 

Remark 2.3 The iteration procedure for two subdomains (2.19)-(2.24) has been introduced in 
[15] for Chebyshev collocation methods, and subsequently applied in [21] to finite element 
approximations. The convergence analysis was carried out in [15] for the Helmholtz operator using 
a separation of variable argument. In section 4, we apply a proof technique based on a suitable 
extension theorem, similar to the one presented in [21] for finite elements. The higher generality of 
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this technique allows us to deal with the case of several subdomains.O 


3 VARIATIONAL FORM OF THE COLLOCATION MULTIDOMAIN PROBLEM 

In this section we give a variational formulation of the collocation multidomain problem 
(2.5M2.9), and, accordingly, of the iterative procedure (2.1 1)-(2.18). The new formulation is 
more suited for the convergence analysis that will be carried out in section 4. 

We confine ourselves to the case of Legendre collocation points, where we will make use of the 
following notation. For functions continuous in we define a discrete inner product, 
approximating that of L 2 (f2 ; ), as follows: 

K N y 

(3.1) (w.z)^ := X X (wz) (x^, y m ) co' k co m , i=l M 

k=0 m=0 


Due to the exactness of the quadrature formula (2.1) for polynomials g 6 P^.j, we have 


(3.2) 


(w,z). N = 


wzdx if wz e P 2N j (Qj) , i=l,...,M 


We also define the discrete bilinear form 


fd w dz \ fd w dz \ 

(3.3) \ N (w,z) := + » % J. N + <«o w An 


3w dz 


For any cp € P° N (T), we denote with 9 e P N (Q) any piecewise polynomial extension of 9 to Q. 

which is determined by the values of <pj j and cp ; solely. 9 can be, e.g., either the discrete 

harmonic extension, or the interpolant extension of 9 . We can now prove the following 
equivalence result. 

Proposition 3.1 The collocation multidomain problem (2.5) - (2.9) is equivalent to look for u N 
such, that := u N1Qi e Pj^Qi), i=i M and satisfies 

(3.4) a u>J (u i>N ,v) = (f,v) i>N Vv€P° N (Q l ) , i=l,...,M , 

u i>N = 0 on SQj n dQ. , 


(3.5) 


i=l,...,M , 
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(3.6) 

u uN “ U i+1.N 

on r. , i=l,...,M-l , 


M 

M M 

(3.7) 

X ( u LN * y) m 

•Ia iN (u jiN ,<?) + I(f 1 ip) k<N V<p€P° N (T) 

j«l k=l 

i odd 


1*1 

i even 


Proof This proof generalizes the proof of an analogous equivalence statement given in [13] for a 
decomposition by two subdomains. 

We start proving the following discrete integration by parts formula: 


Nv 


(3.8) 


a* (w.z) = (Lw,z) iiN + *)(*i • y m ) * (g-jf z ) < x i-i • y«)] ®m 

for i=l,...,M , V w, z e P^f^) such that w=z=0 on n 8Q . 


Indeed, by (3.3) and (3.1) we obtain, using the property (3.2): 


(w,z) = 2 


N y 

m*0 


(sr I) <w-> * + s (y I) <4.y) <*y + <“o », A„ 


>4 




*i-i 


-l 


N y 

— 

m=0 


^a 2 w ' 

y^ z ) 


( x -y m > dx - (f- z )( x i-y m ) + (g^ z ) ( x i-i • yJ 


*i-l 


* 


i 


-X< 

k*0 


-1 


^ 8 2 w n 

va7 y 


(** y) d y + ( a 0 w ’ z >ij< 


Now we get (3.8) using again (3.2). Let now u N be the solution to (2.5)-(2.9). Equations (2.6) 
give: 


(3.9) ( ^ u ij4 » v \n = » v )i jsj V v 6 P°^j (£2^) , i— 1,...,M 

This yields (3.4) by virtue of (3.8). 

We verify now (3.7). For this, let i be an odd integer between 1 and M, and let cp be any function 
of P° N (T).In one has (p = V* 1)+ v/ 0 , where is a polynomial extension of the function 
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r* r* (j) 

which coincides with 9 on C * and vanishes on C l , 1 * j. In particular, we can set y = 0 on Cj" 1 

if 1 * j, j + 1. Using (3.8) we have 


<«) _ 0 ), V r 3u iN 1 

a i^( u iN.v ) = (Lu lN ,v ) iN + 2 -col 9 i (y m ) 

m*0 L lx=Xj J 

V r • 1 ' CO 

2 , 0 ) 10 ) 0 . 11 ^- 01 ^ +-^ (y m ) % (y m ) + (f. ¥ )ijs 

f + ^i+l,N 1 (i) 

- 2, 0) J ca (Lu i+1>N - % - (y m ) <Pi (y m ) + (f.v 

m=0 L ° lx=xj 

(0 r 3u 1+1N 1 o) 

- (Lu i+1N , v ) i+1>N + 2,0)1 9; (y m ) + (f>v )i.N 

m=0 L lx=Xj J 


= (using (2.6)) 

= (using (2.9)) 

) i<N = (using (2.6)) 
(») 

+ < f -Y )i+i,N • 


Then, using (3.8) in £2j +1 we find 

( u uN • V°> - -a i+1>N (“i+l.N • V°) + (f , V 0 )^ + (f , • 


In a s imilar way we find that: 

a U 4 < u ijt • '/ i ' 1) ) = ~ a i-lJSf • v i_1) ) + (f. + ( f . V U 1 ) )i.i.N • 

By summation of the last two equations it follows: 

a iX Kn’ 9 ) = _a i-i (“i -1 V i4> ) (U i+KN , V i} ) 

(3.10) 

+ (f,9) u ,+(fy i * i \ 1J ,+(fy i> ) i+1 jt 

It is now an easy task to see that (3.7) can be obtained by summing the relations (3.10) on all odd i 
between 1 and M. 

We shall now prove the converse, namely that the solution of (3.4)-(3.7) is also solution of 
(2.5)-(2.9). To this end, for each i=l,...,M let P denote a point of Cj. The discrete characteristic 

function associated with P is the polynomial of P N (12 i ) which equals 1 at P and vanishes at the 

remaining points of C ; . Equations (2.6) are obtained from (3.4) taking 9 to be either of the discrete 
characteristic functions associated with either point of Cj*" 1 . Similarly, (2.9) follows from (3.7) 
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taking as 9 a function of P Ny (T.) which equals 1 at a point of C^ 1 , and vanishes at all the other 

points of C*\ then working as done for proving (3.10). Clearly, the relations (2.7) and (2.8) 
follow from (3.5) and (3.6).0 

In view of the Proposition 3.1, it is now clear that the iterative procedure (2. 1 1)-(2. 18) admits 
the following equivalent variational formulation: 


For all j odd ( ) solve 


(3.11) 

a uN( un iN- v )*( f - v )ijt 

V v e P° N (Q ; ) , i=l,...,M , 

(3.12) 

u\ N = 0 

onSQnSQj , i=l,...,M , 

(3.13) 

„ _ f 0 „-l “i-LN + 0 * e „.i) U X 

<> -a.) >41 

on r j , i=3,...,M-l 
onH , i=l,...,M-2 , 


then for all j even. (i=2,...M-l) solve 

(3.14) 

a uN < u \n • v > ■ (f. v) w 

V v € P» N (Q ; ) , 

(3.15) 

u\n = 0 

on d £2 n 8 , 

(3.16) a^u"^ 

.9) ■ - a i-i^( un i-i.N -$) + ( f .9)ij^ + ( f -9)j 

V< PS P’n(T): <P| rj s0 if 


(3-17) a iJv( (u^ , cp) = -a i+1N (u i + i,n > 9) + (f» ^)i.N + ^ 9 e P n (1~) ■' 9| r - 0 if j 5 ** 

For each i odd, (3. 1 1)-(3.13) is a Dirichlet problem, while, for i even, (3.14)-(3.17) is a mixed 

type problem. The above variational formulation of the multidomain collocation iterations 

(2.11) -(2.18) will be used in the next section to investigate the convergence properties 
of this method. 

Remark 3.1 If Q is a plurirectangle, partitioned into subrectangles Qj without internal cross 

points (as in fig.O.l.a), the above multidomain problem and the relative iteration procedure 

(2.1 1) -(2.18) can still be formulated. In such case each vertex of the decomposition lies on the 
boundary of Q, hence the value of each iterate u\n is set to zero (the prescribed boundary data) 
there. 0 

Remark 3.2 For one dimensional multidomain methods using the Chebyshev collocation points an 

equivalent variational formulation is still available (see [14]). However, the bilinear form a iJsf (-,-) 

that holds such an equivalence is much more involved than the Legendre one (3.3). (Actually it 

must take into account the Chebyshev weight functions which blow up at each subdo- 
main interface). The generalization of this equivalence to two dimensional problems is 
not yet available. Thus the convergence analysis we carry out in the next section will 
deal with Legendre collocation points only.<> 
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4 CONVERGENCE ANALYSIS 

In this section we confine ourselves to the case of the Legendre collocation. We prove that the 

sequence (u n i-N ), solution to (2.11)-(2.18) (or, equivalently to (3.1 1)-(3.17)) converges, as n -> 
<», to Ujjj and that the rate of convergence is independent of N. 

For any ye P° N (D we recall that 9 on I~. is denoted by <p i , i=l,...,M-l; then we define 
Rjjj <pe P N (Qj) through the equations: 


(4.1) 

a iJ4 (R|J4 9> v ) - 0 

Vv€ P° N (£2 i ), 

(4-2) 

Rjj,9 = 9. 

onT, 

(4.3) 

Rm9 = 9i.i 

on r, , 

(4.4) 

R uN 9 = 0 

on d £2 0 3 


If i=l or i=M the conditions (4.3) and (4.2) respectivey must be dropped. In view of the analogy 

between this problem and (1.4), <p is referred as the discrete harmonic extension of 9 to Q r 
We now define the following norms. 


(4.5) 


(4-6) 


H<PW= 


• ^ t a uN (^i.N 9 ’ ^i.N 9)1 | 
'-1 

odd J 

f M 


II 9 II 


[ a yj (R^ 9 , Ryj 9)] 


1/2 


> 


i even 


In the next lemma we state that the norms (4.5) and (4.6) are uniformly equivalent on P° N (T). This 

is an extension type theorem for spectral collocation approximations similar to the one stated in [3] 
and [21] for finite elements. Its proof is essentialy based on a result due to Bemardi and Maday 
[ 2 ]. 

Theorem 4.1 There exist two positive constants q and t independent of N such that 


(4.7) 


»9'&d* oll 9l|L “d 11 9 * * 11 9 "ltd • for any 9 e P° (T) . 
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adih 


4g j |o»a 0 p(M) , t»x 0 p(M) , and p (M) * (1 + H" 1 ) max (1, H -1 ) . 

(H ■ min {meas (£2,) , l£i£M) 

Proof Let us fix i such that 2 £ i £ M-l. For any <p e P° N 01 let 9 e P N (£2,) be the solution 
to the collocation problem: 


-AQ UN q> = 0 
Qijt 9 * 9 i 

Qux 9 “ 9 i4 

Aj*9 =0 


inCf“ , 
onT. , 
onT, , 
on 3 £2 n 8 £ 2 ; . 


It is proven in [2] that (for £2, = (-1, l) 2 ) 


1/2 

< 4 ' 9 > 


For a general rectangle £2 it it is easy to see that Cj is proportional to l + meas (£2,) 
Setting v = R^ N 9 - Q uN q> in (4.1) we obtain 


(4.10) 0^ij4 9 * Qij< 9 » ^i,N 9 • Qij4 9) - ‘^n 9> ^ij4 9 * 9) • 

On the other hand we Tecall that (see [12] and [ 6 ,Ch.ll]) 

(4.11) q II w II 2 , £ a, N (w,w) V w € P N (£2,) 

H (tij) u " * 

and 

(4.12) I a LN (w,z) I S C, II w II , II z II . Vw,zeP N (£2) 

where Cj and C 3 are positive constants independent of N. Hence, from (4.9) and (4.10) we 
deduce 

<413 > M- <*» * • *W « s C« [• 9, ^ - « 9,, ^ ] 
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In order to estimate the right hand side of (4.13) we use now the continuity of the trace operator 
from H‘ (Qj) to H ,/J 0 fl,) (see [18]) 

(Ryt 9 . *iN 9) * Cs [ n ^tjt 9 , H , ( Q i+ , ) + 11 * ^tpj. j} ] 

where C5 is proportional to p (M) 

Using again (4.1 1) we also get 

(4.14) [a^CR^ 9, 9)] 2 S C 6 (la^ 1(N (R j+l<N 9. R i«.i.N 9)] 2 + 9. R j.i,N 9)1 2 ) • 


For i=l or i=M working in the same way we have 


(4.15) 


l a l.N ^l.N 9> ^l.N 9 )] 2 ^ C 7 la^ (Rjjj 9* 9 )] 2 » 

l a M.N 9» ^M.N 9 )] 2 * ^8 f a M l.N ^M-l.N 9 * R Mi.N 9)f • 


Finally summing up (4.14) and (4.15) for i odd and for i even we deduce the desired result (4.7).0 
In order to prove that the solution u"^ of (3.1 1)-(3.18) converges, as n -* •», to the one 
of (3.4)-(3.7), we look for the operator T e such that 

(4.16) * u lN * T e (u^ - u^) on f , i even . 


The operator T e can be defined as follows. For any <p € P° N (T) let w^. € P N (£2,). i=l,...,M, i 
even , be defined through: 


(4.17) 

a i.s ( w ln * v > - 0 

Vv€P e K (Q,) 


M 

M 

(4.18) 

2* a^ (Wj* , R i>N v) = - 

i« j 

X a uN (R^ 9 . R^n V) V v e P° N O') 

is] 


i even 

i odd 

(4.19) 

< 

H 

O 

on 3 O r> 9 Clj 


Then we define 

(4.20) Ttpr-wyj for i=l,...,M-l , 

and 

(4.21) 


T e 9 :* 9 Ttp + (1 - G) 9 . V96 P“jj (T) . 
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Remark 4.1 Due to (4.1) and (4.17) one could substitute in (4.18) with any other extension 
of V in £i ; . without affecting the definition of w^.O 

Remark 4.2 The operator T can be interpreted as follows (suppose M > 2). For a given <p € P° N 
(T) and an even number i with 2 £ i £ M-l one can find cp and R i+1 N <p by solving the 
Dirichlet problems (4.1)-(4.4). Then compute the left flux at the boundary side F x and the right 
flux and the boundary side F as follows 

Di-ljn = £ R i-l ji <p < X N » yj + LR i-l <P ( X N - y m) “ N 

D i+ljn = ' ^ ^+1 Jt ( x 0 +1 > y m ) + LR i+l *P <*0 1 ’ ym) “o 

for m=l,..„ N-l. The functions w^ are now obtained by solving the mixed problem 

’^N € p n ( Q i) 

Lw^ n = 0 in Cf 

dw i S ■ i 

. ( X N> y m ) + Lw iJJ < X N’ ym) = D i+l,m for m=1 ’ - ,N-1 

3Wj M i 

-g — y m ) - Lw UN ( x 0- y m ) <°0 = D i-ljn fOT m= l* - >N -1 

w iJ4 = 0 oncf . 


Finally Tcp on F j and F is obtained by computing w^ on F } and F. 

The operator which associates to the values of Dj. ljn and D i+1 ^ the values of w^ N on Fj and F 
is a discrete Poincar6-Steklov operator (see [1]). 




r 

r 




i-l 

i 


r 

R <p 

w 

R © 

r 

i-2 

i-l.N 

i.N 

i+UN 

i+i 


no n 

i-l i i+1 


FIG.4.1 Computation of T <p 
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In the following theorem we prove that for suitable values of 0, T 0 is a contraction. This, together 
with (4.16), will imply the convergence of u"^ to u^ N for n -» <». 

Hlggrsm.4,2 Eor_each M S 2 there exists 9* > 0 such that for anv N > 0 the following holds 

(4.22) V 9 € (0, e*) 3 k (0) < 1 : II T 0 <p 11^ S k (0) II <p 11^ V <p e P° N (T) 

Moreover there exist 9’. 0" and k mill 0 < 0' < 0” < 0* and k < 1 such that for all N > 0 

(4.23) V0€ (0', 0"] k(0)*K<l . 

Proof From (4.5) and (4.21) we have, using the symmetry of the forms (•, ) 

2 M 

(4.24) II T 0 <p ll^j * 0 II T<p II £* + 20 (1-0) X a^ (R^ cp, R. N T<p) + (1-0) 2 II cp ll 2 ^ . 

i*l 
i odd 

Moreover, using (4.17) and (4.18), we have 

M M M 

(4.25) X a lN (R ifN <p, R i>N Tcp) = - X a i>N (w^ , R Tcp) = - X (R. N T<p, R^ T<p) , 

1*1 i*l isl 

i odd i even i even 

hence, from theorem 4.1 and (4.24) we deduce 

(4.26) II T 0 <p ll 2 ** < 0 2 o II T(p ll 2 even + (1 - 0) 2 II tp ll 2 ^ - 20 (1 - 0) H Ttp ll 2 even . 

From (4.25), the Cauchy-Schwarz inequality and (4.7) we also have 

(4.27) IIT<p<Us/7ll<pl U . 

In order to find a lower bound for II Ttp ll even we remark that from (4.17)-(4.18) we have the 
fundamental relation 
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M M M 

X (Ry< 9* ^tN 9) ■ * X a tjj ( w ijj • ^ijj 9) = * X ajjq (R^ T«p, R uN <p) , 

i»l i*l i*l 

i odd i even i even 


Using the Cauchy-Schwarz inequality, the property (4.7) and the definitions (4.5) and (4.6) we 
deduce: 


(4*28) H 9 "odd ^ 11 T 9 "even • 

Jx 


From (4.26),(4.27),(4.28) it follows 


, f 2 2 , 20 (1-0)) 2 

'V&d* 0 o +0-8) »9l&d • 


) 


If we define 


k(0) := 


2 2 

0 (a x + x + 2) - 20 (x + 1) + x 


-.1/2 


We can readily see that (4.22) holds and that 

k (0) < 1 iff Ck0<0* = min 


2 (x+ 1) 

1 __ 

ax + x + 2j 


Then (4.23) follows from the continuity of k(0).O 


We can now derive the following convergence result. 

Corollary 4.1 L& u i>N , i»l,..., M, be the solution of (3.4)-(3.7) and let u n iiN , for n > 1 aM 

i«l,... , M, be the solution of (3. 1 1)-(3. 18). If 8„e [0’, 0"] we have 

(4.29) 

moreover 

(4.30) 

where C is a constant dependent on e° but independent of both N and M. 

Proof The estimate (4.29) follows using repeatidly (4.16) and theorem 4.1. Due to the theorem 
4.1 and to (4.11), the estimate (4.30) can be deduced from (4.29).0 


n+1 


il (u^ 1 - Ujjj) ir Hodd S K II (u°j^ - Ujj^j. IIq^ 


M 


,1/2 




SCic 


n+1 


Remark 4,3 In view of (4.8), if H is much smaller than one, then a and t grow like H' 2 . If 8 
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T+ 1 


is the value of 0 that minimizes K (0), one has 0^ 


~2 and k (0^) behaves like 

o t + x + 2 


(1 - CH 2 ) 1/2 , where C is a positive constant 0 

Remark 4.4 The proof of theorem 3.1 and corollary 3.1 can also be applied to the case of a 
plurirectangle £2 partitioned into rectangles £2j without internal cross points (see fig.0. 1 and remark 


3.1).0 


5 MATRIX INTERPRETATION OF THE ITERATIQN-BY-SUB DOMAIN METHOD 

To solve the multidomain problem (2.5)-(2.9) we have proposed, in section 2, an iterative 
method based on a sequence of differential problems to be solved in each subdomain. An (a priori) 
different point of view consists of applying the influence (or capacitance) matrix method for 
solving (2.5)-(2.9). This approach, in the framework of spectral approximations, has been 
extensively pursued by Peyret and his Coworkers (see e.g. [17], [24]) and by other authors (see 

[20], [22]). The influence matrix coincides with the Schur complement of the matrix of 
the system (2.5)-(2.9) with respect to the interface variables (see e.g. [9]). It is precisely 
the matrix of the system of the interface unknowns and it is derived from the global 
system by block Gaussian elimination (see e.g. [3]). In this section we show that the 
iteration-by-subdomain method introduced in section 2 is equivalent to a preconditioned 
Richardson iterative method for the resolution of the influence system. In particular we 
will give the precise structure of the preconditioner for the influence matrix. 

Wc proceed now to construct the influence matrix associated with the problem (2.5M2.9). 

Let [Tjj 6 P° N (T.), i =1,.„,M-1] be M-l given polynomials. Correspondingly, for each 
i=2,...,M, denote by v^ the solution to the problem 


6 (®i) 

Lv in = 0 

in Cf“ 

FT* 

II 

in C^' 1 

< 

£ 

II 

O 

inC Fi 

O 

II 

A 

>~ 

in cf 


Moreover, for each i=l,...,M-l we let be such that 
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0 0 
FIG.5. 1 Boundary conditions for and 


Let us denote with ri- the vector of the values attained by the polynomial T). at the collocation points 
CF of r. Then we define S'j, S' 2 , S‘ 3 , S ' 4 to be the square matrices of order N-l such that 

(5.3) (S', q.) m := ^ (x' N , y J + Lz^ (x' N , y m ) co. 

(5.4) (S‘Tl, l)m := v^ (x'o , y m ) + Lv lN (x' 0 , y m ) ( 0 ^ 

(5.5) (SJ T|.) m := ZuN ( X q , y m ) + Lz i>N (x‘ 0 , y m ) co^ 

(5.6) (Si Ti M ) m := ^ v^ (x' N , y m ) + Lv iJ( (x‘ N , y m ) C 0 j 

for m=l,...,N-l. S'] is defined for i=l M-l, S ’ 2 for i=2,...,M and S 3 , S ' 4 for i=2,...,M-l. 

We shall call for brevity interface flux on T. of a function w defined on the quantity 

^ w ( x n « y m ) + Lw ( x n • y m ) ®i for m=1 N - 1 • 

Then the matrix S*j represents the operator which associates to r). on T. the interface flux on T. of 
the left discrete harmonic extension of rj. on S' 2 , S 3 , S ' 4 can be interpreted analogously. 


I 



FIG.5.2 Interpretation of S l j, S* 2 , S^, S l 4 . 


We now define the block tridiagonal matrix S, of order (M-l) (N-l), as follows. 


2 2 3 

s s + s 

4 1 2 


o 


o 


M-2 M-I 

M-l 

s + s 

s 

1 2 

3 

. M-l 

M-l M 

s 

s + s 

4 

1 2 


We claim that the matrix S in (5.7) is actually the influence (or capacitance) matrix of the system 
(2.5)-(2.9). As a matter of fact, let us solve the homogeneous collocation problems 


\n 6 p n < Q i) 

(5.8) | Lu; - f inCf 


.u^ = 0 in C ‘uC b 
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for i=l,.— »M and compute the vector b * [b 1 ,...,b i ,...,b M .i] T defined by 


(5.9) (b,) m = |;\n (x‘ N , y m ) + (Lu^ - f) (x‘ N , y m ) ©. - ^ u i+1N (x” 1 , y m ) + 


(Lu i+1J< * 0 (*o +1 1 • y«n) ®i * for m=l,...,N-l, 


Then define the vector % = _j] T as the solution of the linear system 


(5.10) 


s 5 =b 


It is now readily seen (using (5.1)-(5.10)) that the solution of (2.5)-(2.9) is given by 


(5.11) 


'«U< = “iJt + w uN + z uN 
• u ljt = u l Jf + Z 1 .N 
• u MJt 31 u MJt + W MJ< 


for i=2,...,M-l 


where w^ and z i<N are obtained from (5.1) and (5.2) taking rj- = Thus, the matrix S defined 
in (5.7) is actually the influence matrix associated with the problem (2.5)-(2.9). 

Remark 5.1 It can be proved, using repeatedly the integration by pans, that in the case of the 
Legendre collocation the matrix S is symmetric. On the other hand S is not symmetric if the 

Qiebyshev collocation points are used. S is again non symmetric in the case of both Legendre or 

9Ujj^ ^i+l,N 

Chebyshev collocation if the pure flux condition ^ — on T. is used instead of (2.9).0 

The relations (5.8)-(5.11) suggest a possible algorithm, based on the influence matrix 
system, for the solution of the multidomain problem (2.5)-(2.9). The crucial step of this algorithm 
is the use of an efficient method for the solution of the system (5.10). To this end, preconditioned 
iterative methods (which do not require the explicit construction of the matrix S but simply the 
calculation of the product of S time a given vector) are currently used (see e.g. [3], [4], [5], (7], 

[8], [11], [16]). In this section we show that our iteration by subdomain method amounts 
to the use of the Richardson iterative method for (5.10) preconditioned with a proper 
matrix P (see (5.15) below). 

Let us recall that if P is any non singular matrix of order (M-l) (N-l), the Richardson iterative 
method for solving (5.10) preconditioned with P is: 

(5.12) P ( £ n+1 - 1") = 0 n (b - S|") . 

The iteration matrix at the step n is 

(5.13) (1 - e n ) I + 0 n (I - F 1 S) . 
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Let us define the matrix P as follows (we assume, e.g., M odd). 



-2 2 

_ 


S S 

2 3 

s 2 s 2 

4 1 

s 4 s 4 

2 3 

4 4 

o 

(5.14) P = 

s s 

4 1 

M-l M-l 

s s 

2 3 


o 

M-l M-l 

s s 

4 1 —1 


We denote with T e and T the matrices associated with the operators defined in (4.20) and (4.21) 
and bearing the same name. 

proposition 5.1 The following identity holds 

(5.15) T e »(l -ejI + ea-P * 1 S) . 

Proof In view of (4.20) and (4.21) it is sufficient to prove that T * I - P* 1 S or equivalently that 
PT » P-S. Recalling (5.3M5.6) the action of P consists of solving a sequence of Dirichlet 

problems in ft. for i even, and then to compute the interface fluxes on T. Recalling the 
interpretation of the operator T given in the remark 4.2, it is readily seen that applying the matrix 
PT to the values of a polynomial on amounts to solve a sequence of Dirichlet problems in ftj 

for i odd with such boundary conditions, and then to compute the interface fluxes on T of the 
obtained solution. 

On the other hand the matrix P-S is given by 
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l 


P-S = 





-S 


M 

2 


From (5.3)-(5.6) we deduce that the matrix P-S is associated with the solution of Dirichlet 

problems in Qj for i odd and to the computation of the interface fluxes on T. This shows that 
PT=P-S. 0 

Remark 5.2 (The case of two subdomains') If M is even it is possible to give an analogous 
interpretation of the iteration-by-subdomain method. In the particular case in which M=2 we have 

(5.16) S = S 1 i + S 2 2 and P = S 2 2 , 


and therefore (5.15) becomes T e = (1 - 0) I - 0 (S 2 ^)' 1 S^). If fij and Qjhave the same measure 

and L = -A + al, with a e R + , then S 2 2 = S* x , hence taking 0 n = 1/2, Vn, yields exact 
convergence in two iterations (this result was already found in [15]). 0 


Remark 5.3 (Condition number of the preconditioned matrix) Let us consider the case of the 
Legendre collocation with M=2. Recalling the definitions of S 1 j and S 2 2 (see (5.3) (5.4)) and 
using the discrete integration by parts formula (3.8) we deduce from (5.16) that 

(5.17) Vcp,f e P^ (17) (S <p, y) Nr = a 1>N (R lN <p> Ri,n V) + ^.n ^2,n *P* ^2,n V) • 
where <p is the vector such that <p m = <p(y m ) m=0,...,N y , and 


(5.18) 


Ny Ny 


< S M) NXi II^9 k y m o> m 

l m=Uk=0 


! 


Similarly the preconditioner P = S 2 2 satisfies 



(5.19) 
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V <P» v € Pn (Tj) (P9* V) N rj = a 2jvf (®2JJ V* ^ZN v) • 

From (S.18),(S.19) and theorem 4.1 we obtain the fundamental relation 
(5.20) V <p e Pj (Tj) (P 9, 9) NJ] ^ (Sep, 9) n >r 2 1 (P9, <P ) N x 


where t is a constant independent of N. This ensures that the eigenvalues of the matrix P**S are 
real and positive. Moreover the ratio X maY A rn j n between the maximum and the minimum 
eigenvalue of P**S is bounded by x.0 


Remark 5.4 (Choice of the relaxation parameter 8 n ). The interpretation of the iterative 
method(2.U)-(2.18) as a Richardson scheme suggests some convenient choices of the parameter 

6 n appearing in (2.18). When information is available on the extreme eigenvalues of the precondi- 
tioned influence matrix the optimal parameter 0 = can be used for all n (see e.g.[7]). 

X + X . 

mix min 

Otherwise a dynamical choice inspired by a "minimal residual” strategy can be used (see [1] for the 
case of multidomain finite element approximation). 

Another strategy is suggested by the proof theorem 4.2. The minimal value of the contraction con- 

2 

stant k(0) is achieved for 0^ * (x + 1/c x + x + 2). In [21] the values of 0 R are choosen so that 

sup {0 n J = 0^ and no extra computational work is needed to obtain 6„. The same strategy 
can be applied to the case of spectral collocation approximations here discussed. <> 

6 OTHER PRECONDTTIONERS 


Based on the interpretation of the preconditioned Richardson method for solving the 
influence system (5.10), one can formulate other iteration-by-subdomain procedures by simply 
taking a preconditioner P in (5.12) different than the one in (5.14). For instance if we take 


1— - 2 


(6.1) P = 


o 



S 


M 


2— J 
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then the iteration (5.12) yields the following iteration by subdomain procedure. (We identify, for 
simplicity of notation, the vector with the polynomial e p° N (T) whose values at the 
collocation points are the components of £ n ). 

Given £°, look for U^e P N (fty for i=l,...,M-l such that for any n £ 1 


( 6 . 2 ) 


LU* » f 


m<^“ , 
n-1 r ; 


U? -§j‘ inC ' , U?=0 inCfand,ifi£2 , 


n-1 


«r -v. fee 


r : 


i-l 


Then, for i=2,...,M look for V"; 6 P N (Q ; ) such that 


(6.3) 


LV?*f inCf“ 

Jjj V* - (LV“ - 0 » , - 4 U V d-ur, • 0 » C r ‘-' 

lvf.< inC Fi , Vf = 0 in c|* if i < M, and = 0 in 


Finally is defined by 

(6.4) £ « 0 n V n M + (1 - 0 n ) on C l for i=l,...,M- 1 

Once I;”' 1 is available, one can compute IPj, i=2,...,M by solving the M-l independent problems 
(6.2). Then the M-l independent systems (6.3) can be solved simultaneously. At each step of this 
algorithm one has to solve 2(M-1) differential problems, instead of M as prescribed by the method 
(2.1 1)-(2.18), however in this case the degree of parallelism of the algorithm is M-l, instead of 
M/2. 

The choice of the preconditioner 
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(6.5) P = 



gives rise to an algorithm analogous to (6.2)-(6.4). In this case the Dirichlet problems are solved in 

Q} for i=2,...,M, and the mixed Dirichlet-Neumann problems in Cl ; for 

Another preconditioner which corresponds to an iteration-by-subdomain technique is 


S 

S 


l 

1 

2 
4 



1 


(6.6) P = 





S 


rj 


The differential interpretation of the preconditioned Richardson iteration is now the following. 

First one solves M-l independent Dirichlet problems in i=2,...,M, with £ i n ' 1 as boundary data. 

Then computes the interface flux <J>. j on ^ of the solution in £2 ; . Finally one solves sequentially 

M-l mixed problems going from £2j to Q M . r For the mixed problem in fl, one enforces the flux <{>. 

in r., while on I", j the value of the solution coming from Qj_j is enforced as Dirichlet data. We 
note that these mixed problems are not solvable in parallel. 

An analogous algorithm is obtained by using the upper triangular preconditioner 
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The convergence analysis for these preconditioners is under investigation. 

Notice that in the iteration by subdomain methods we have considered so far the subdomains can 
be handled in parallel whenever the preconditioner P has a block diagonal structure. The degree of 
parallelism is precisely the number of diagonal blocks. 

7 DECOMPOSITION OF Q WITH INTER NAL CROSS POINTS 

In this section we consider the case of decompositions of ft with internal cross points (see 
fig.0.1 b). 

For simplicity we consider the case of a square ft partitioned into four squares ft^ i=l 4 (see 

fig.6.1). The general case can be studied combining the results of this case with those of the 
previous sections. 



F1G.7.1 The decomposition of ft. 
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4 


Wc denote by <X> the space of the functions which are the traces on F := vj r. of the functions of 

j=l 

HVfi). The variational multidomain formulation of problem (1.1) is given by (1.5)-(1.8) 
provided Uj is identified with u 5 . 

The collocation points in , i=l,..., 4, are defined in section 2. We recall that here C*~‘ is the set 

of the internal collocation points of the interval T., i=l,...,4. The cross point 0 is the unique 

collocation point common to all the subdomains. For simplicity we look for a spectral solution 
which is of equal degree (N) in the four subdomains. Hence the weights defined in (2.2) and (2.4) 

are such .hat o>‘ k =to k ,Vk = 0,...,N, moreover we have co^ = co Q , and we denote this value with 
to. 

Let us consider the following multidomain spectral approximation of problem (1.5)-(1.8). For 


i=l,... 

,4 we look for u^ e P N (Q;) such that for i=l, ... ,4 (we identify u J>N 

with u 1J4 ) 

(7.1) 

Lujjj = f in Cf 


(72) 

r £ 

u uN = u i+i ,N “ C 


(7.3) 

Uln = 0 hi 


(7.4) 

9u i,N 3u i + l.N „ A _ . . Tj 

= ‘(Lujv - f)co - (Lu i+1J< - f)oi ui C 


(7.5) 

A 

r 9u 4J4 5u 4.N ) 

t“5r +_ 5rJ (0) 


4 

- - X (Lu^ - f) (0) to . 

i*l 



The u’s are the outward normal directions to r ; , ordered clockwise as indicated in fig.6.1. 

The conditions (7.4) are analogous to the ones in (2.9). The "cross point condition" (7.5) gives a 
relation between the four polynomials u ijs( , i=l, - ,4. In the spirit of the method (2. 1 1 )-(2. 1 8) it is 
possible to formulate an iterative method for solving (7.1)-(7.5). The idea is to solve, at each 
iteration, two Dirichlet problems in Cij and £2 3 , and two mixed problems in flj and £2 4 as follows. 

For i=l,...,4 let g® e P N (Tp be such that g® = 0 on and g°i = g °2 = g °3 = g °4 in 0; we look 


for a sequence u \n> n2 1, satisfying: u”^ e P N (Cl), i=l,...,4 and 

iiifljipd-fly. 

(7.6) LuJj, = f in (**, = g^ 1 in C l , u?_ N = g^ 1 in C 4 , u" N = g?' 1 in 0, u" N = 0 in C* 
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(7.7) Lu^ » f in Cf* , u£ N = g^ 1 in C 2 , u^ = g£' 1 in c \ N = g" 1 in 0 , u"^ =0 in C$ 

iiL&^aDsLC^: 

(7.8) Lu£jj = ^ *® p2* * U 2J< = ® *° ^“2 

(7.9) . —H: + (Lu^ - f)<0 = - — — - (Luu< - f)c» inC 


(7.10) 


*>! 

<h>„ 


du. 


- (Lu^j - f)o + (Lu 3 - f)co in C 


dv. 


(7.11) Lu 4 jj = f in Cf , u 4 ^ = 0 in 


3 uJ u u r 

( 7 . 12 ) — --(Lu 4 jr f)a> = -— + 0-113^-00) inC 

<h > 3 <n > 3 


9iij u Ou, m r. 

(7.13) — — + (LU 4 J 1 - f)tD = - — (Lu 1<n - f)<o inC 

5u 4 dx> 4 

(7.14) - (0) + + ^r](°) + ( Lu In * 0 (°>“ + - 0 «»“ = 

= - (^r - ) (0) + (tgt * TF ] (0) - ^ - f)(0)to - (L u 3.n - o (O)co 


(7.15) u^ = inO . 

The polynomial g"} are defined recursively as follows: 


(7.16) 


8? - e , “In ♦ ■ 9 „> 8?' 1 

inC** 1 

4-8. <■«*«-»,) 8?' 

inC* 

«S * \ u !.n + (‘ • e „) e!’ 1 

mC 3 

<n • 0 „) C 

mC‘ 


for i=l 


Note that problems in Qj and are independent, while those in and Q 4 are coupled through 
the two cross point conditions (7.14) and (7.15). 

If the Legendre collocation points are used, the problem (7.1)-(7.5) can be written in a variational 
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fonn which coincide with (3.4)-(3.7) provided M=4, u J>N is identified with u,^ , and P° N (I~) is 
substituted with the space 

4 

<I> N ;* ( 9 : u rj -» R : 9 6 P N (T.) , i=l,...,4 , 9 is continuous in 0, 9=0 in 9C1 n (U T.)}. 
i -1 1 

We skip the proof which is analogous to the one of proposition 3.1. With this variational 

formulation we can now follow the guideline of section 4 to prove the convergence, as n —♦«*>, of 
the iterative scheme (7.6)-(7.16). In this case however the extension theorem 4.1 takes a different 

form. Precisely (4.7) holds now with two constants o and t which depend on N as follows 

o s Oq log N, x = XpIogN . 

A consequence of this is that the constant 6 * of theorem 4.2, as well as the error reduction factor 

k(0), for 0 € (0, 0*), depend now on N. The convergence result (4.29) is still true, but k tends to 

one as N — > 00 with a logarithmical speed. The precise form of this result and a detailed proof of it 
can be found in [25]. 
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